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1 .  INTRODUCTION 


1 . 1  Purpose  of  the  Study 

This  study  was  based  on  the  hope  that  we  could  utilize  the  power  of  the 
computer  to  both  design  and  construct  pattern  recognition  holograms.  At  the 
beginning  of  this  study  the  only  computer  designs  were  limited  to  80  picture 
elements  (pixels)  in  the  hologram  plane1  and  computer  hologram  writers  were 
limited  to  a  few  million  total  pixels  with  positional  accuracy  of  roughly  one 
part  in  a  hundred  (or  worse)  in  each  direction.  Accordingly,  our  goals  were 
to : 

o  Computer  design  holograms  of  arbitrary  size  (number  of  pixels) 
optimized  for  any  selectable  goal  and 

o  Construct  a  computer  hologram  writer  capable  of  one  part  in  a 

million  location  accuracy  and  capable  of  writing  108  to  109  pixels. 

1 . 2  Results  of  the  Study 

Both  goals  were  met.  As  a  result  we  can  now  design  and  construct  the 
optimum  hologram  for  any  reasonable  task.  This  appears  to  be  a  revolut ionary 
advance  over  the  prior  art  in  which  some  ad  hoc  optimization  of  optical 
holograms  was  possible, ^  but  no  systematic  approach  to  full  optimization  in 
either  design  or  construction  was  available. 

As  this  work  progressed,  it  became  obvious  that  computer  optimization  of 
optical  masks  was  possible  in  principle  and  that  the  optimization  mathematics 
(eigenvector/eigenvalue  solutions)  was  subject  to  optical  solution.  We 
believe  that  these  two  realizations  are  important  but  not  key  to  our  primary 
results.  Accordingly,  they  are  relegated  to  appendices. 
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2.  HOLOGRAM  DESIGN 


2.1  Introduction 

The  use  of  statistically-derived  optical  filters  has  been  discussed  by 
many  authors  over  the  past  several  years. *-lt  One  of  the  methods  discussed  was 
called  the  "generalized  matched  filter"  or  GMF  because  it  reduces  to  the 
matched  filter  for  simple  enough  situations  (very  simple  invariant  objects  in 
colored  noise)  the  GMF  gives  different  results.  The  basic  idea  is  to  view  the 
M  x  N  complex  feature  vector  and  use  the  classical  statistical  method  of 
linear  discriminant  analysis  to  find  that  linear  combination  of  the  MN  samples 
which  best  separates  (in  an  expected-value  sense)  the  class  of  images  of 
interest  from  all  other  images  and  noise.  Thus  all  "uninteresting"  images  of 

interest  from  all  other  images  and  noise.  Thus  all  "uninteresting"  images  and 

noise  are  lumped  together  according  to  their  probability  of  occurrence  to  form 
the  set  of  images  against  which  we  seek  to  discriminate. 

The  basic  Ideas  are  easy  to  express  mathematically.  Any  Input  scene 
f(x,y)  will  produce  a  Fourier  transform  pattern  F(u,v).  In  discrete  terms, 
the  Fourier  transform  samples  can  be  written 

F  *  F(mAu,  nAy),  (1) 

mn 

where 

0  >  m  >  M, 

0  >  n  >  N, 

Au  •  sampling  spacing  in  the  u  direction,  and 

Av  ■  sampling  spacing  in  the  v  direction. 
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This  reordered  F^  matrix  becomes  the  vector 


G  =  (Gj,  G2,  ....  G^) 


We  then  seek  an  MN  component  vector 


c  ‘  V  S . cm  > 


such  that  the  numbers 


(3) 


(4) 


P  =  C 


C1G1  +  C2G2  + 


+  C  G 
MN  MN 


(5) 


separate  into  one  set  {Pj}  for  the  desired  set  of  targets  and  another  set 
{PT}  and  {PcnI  are  "optimally"  separated.  Our  definition  of  optimality  is 
chosen  for  mathematical  tractibility.  We  choose  C  such  that 


vr 

C _ BC 

2T  w6 


(6) 


is  maximized.  The  "between-class  scatter  matrix"  B  is  given  by 


B  -  I  N.  [5.  -  <6  ]  [2.  -  <&  l1 
i  =  l 


(7) 
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where 


*  number  of  entries  in  { PT } > 

*  number  of  entries  in  {P^}, 

Gj  *  mean  G  for  {P^. } , 

G2  ■  mean  G  for  |pCN}»  and 
<G>  »  mean  G  for  the  entire  data  set. 

The  "within-c lass  scatter  matrix"  is  given  by 


u  ■  W  +  w 

-  1  -2 


(8) 


where 


“i-J,  <8i,i  -  5i><5i.i '  V1 


.  -  iCh/  G  vector  in  {P_ }, 


(9) 


*2 


1,  <02,r  G2KS,i  -  V  ■ 


i-1 


C2  i  -  ith/  £  vector  in  {PCN}  • 


(10) 


Differentiating  Equation  (2-6)  with  respect  to  C  leads  to  a  generalized 
eigenvalue  problem  of  the  form 
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[b  -  aw]  c  =  0  (11) 

The  eigenvector  with  the  highest  eigenvalue  X  is  the  C  whose  components  define 
the  CMF. 

Clearly  a  256  x  256  CMF  requires  solution  of  a  generalized  eigenvalue 
problem  with  65536  x  65536  B  and  W  matrices.  Our  initial  computer  program 
allowed  us  to  solve  80  x  80  problems,  e.g.  and  8  x  10  CMF.  We  know  of  no 
general  methods  to  solve  such  large  eigenvalue  problems.  The  primary  goal  of 
this  paper  is  to  observe  that  the  ratio  of  the  determinant  of  the 
between-c lass  matrix  is  maximized.  Here 

E  [p  ]  =  expected  value  of  P  for  an  input  in  set  (x }  where 
x  =  T  or  CN  and 

o^  =  standard  deviation  of  P  values  for  inputs  of  set  {x'. 

In  words,  we  want  to  optimize  the  T/CN  separation  relative  to  the  tightness  of 
their  clustering.  A  more  precise  statement  is  available  elsewhere.5 
The  vector  C  is  then  retranslated  into  a  two-dimensional  format.  The  mn 
component  of  C  is  the  transmission  of  the  GMF  at  mAu  and  nAy. 

The  purposes  of  this  paper  are  to  show  that: 

(1)  Very  large  MN  values  can  be  handled  quite  easily  if  we  arbitrarily 
assume 


Elpk-  PV  *  <pt-  P»>1 


k* 


(12) 
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where 


E(Pk) 


E(V 


6k , 


fl  if  k  -  l 
>0  if  k  *  t 


and  *  indicates  complex  conjugation; 

(1)  The  calculation  problem  can  be  simplified  tremendously  by  assuming 
that  the  individual  components  of  G  have  no  cross  correlations  and 
hence  that  J5  and  W  are  diagonal, 

(2)  This  assumption  leads  to  very  useful  GMF's  relative  to  some,  not 
necessarily  all,  of  the  alternatives, 

(3)  The  GMF  correlation  peaks  so  calculated  are  one-pixel  wide  (as 
opposed  to  the  auto-correlation  pattern  produced  by  the  MF  or 
matched  filter  which  is  twice  the  diameter  of  the  object). 

Thus  this  is  intended  to  demonstrate  GMF  calculation  and  operation.  A 
fair  and  statistically-useful  comparison  of  all  the  available  pattern 
recognition  filters  would  be  desirable  but  is  well  beyond  the  scope  of  this 
paper . 

2.2  Digital  Generalized  Matched  Filter  Demonstration 

To  demonstrate  the  digital  GMF  we  used  a  binary,  hand-drawn  silhouette  of 
an  airplane  on  a  128  x  128  grid.  Figure  2.1  shows  a  display  of  the  airplane 
on  a  Raratek  graphic  display  system.  The  nature  of  the  Ramtek  display  is  such 
that  a  square  input  is  displayed  as  a  rectangular  output  with  14  units 
horizontally  for  each  10  units  vertically.  We  will  use  this  display  to 
produce  color  and  black  and  white  representation  of  the  images  generated 
during  our  digital  experimentation.  We  have  chosen  to  concentrate  on 
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reduction  of  within-class  variability  of  filter  response  as  a  first 
demonstration  of  two-dimensional  GMF's.  Therefore,  we  worked  with  both  scale 
and  rotation  variations  of  that  input  pattern.  The  pattern  was  magnified  and 
rotated  bv  digital  means  through  a  range  of  magnifications  and  rotation  angles 
indicated  in  Figure  2.2.  Figure  2.2  will  be  important  because  we  will  use 
that  same  display  format  to  compare  CMF  with  with  matched  filter  results  later 
in  this  pape r . 

To  produce  the  generalized  matched  filter  we  use  as  the  set  of  objects 
that  we  wish  to  put  in  the  same  class  the  nine  magnification  and  rotation 

combinations  of  Figure  2.2.  We  created  a  GMF  which  best  separated  this  set  of 

airplanes  fnim  a  second  set  of  inputs  comprising  Gaussian  noise  of  zero  mean 
and  standard  deviation  equal  to  that  of  one  tenth  the  magnitude  of  the  peak  of 
the  Fourier  transform  of  the  unit  magni f icat ion,  zero  rotation  Fourier 
transform.  We  then  generated  the  generalized  matched  filter  using  at  each 
point  in  the  128  x  128  Fourier  transform  plane  the  nine  Fourier  transform 
components  as  well  as  the  128  x  128  array  of  white  noise  samples  as  described 

in  detail  in  the  Appendix.  The  filter  was  generated  by  assuming  that  there 

was  no  cross  correlation  between  any  two  points  in  the  Fourier  transform 
plane.  This  is  demonstrably  not  the  case,  but  handling  the  general  case  of 
covariance  matrices  with  off-diagonal  term  and  the  very  large  sizes 
encountered  in  our  problem  would  be  impossible  with  existing  equipment.  It 
would  involve  solving  an  eigenvector  problem  of  the  form  of  Equation  (6)  for 
the  ( 1 2 8 ) 2  points  in  the  Fourier  transform  plane.  Our  simplification  allows 
us  to  treat  the  problem  in  conveniently-sized  pieces  as  detailed  in  the 
appendix. 

We  shown  in  Figures  2.3  and  2.4  horizontal  and  vertical  traces  through 
the  correlation  points  produced  by  using  the  matched  filter  to  the  M  =  1,  0  = 
0°)  scene.  Figure  2.6  shows  the  magnitude  squared  of  the  matched  filter. 
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Figure  2. 
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Figure  2.8.  Horizontal  Traces  Through  the  GMF  Correlation  Planes 
Corresponding  to  the  Various  Inputs  of  Figure  2.2 
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.9.  Vertical  traces  as  in  Figure  2.8 
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We  ..in  now  .  .  > !  ti  ;  > . » r  -  ■  these  with  the  MF  r--->  Mr  J.-si  gned  to  group 

a !  ’  nine  of  t  ti«  •  images  i  nt  '  one  .'Its-  is  ->!i  vi  horizontal 

and  wrti  .i!  t  r.t.vs  thr’ugh  the  corr.  .  i!  i  :>•  i,  -  !'•  -h  "w:  i  n  Fi  gur.*s  2.8 

and  2  .  ^ .  r  i  i r * •  2  .  ’.**  sh.'ws  1 1 s •  ■  "c^rr--  1  -it  io>;  :  .  r:-"  :>r  •  :  i  •  • :  t  he  n  =  1,  r-  = 

'  i  :i ;>•  1 1  md  titi*  MF. 

It  is  i  :i  t  s  t  i  n  it  t  '  o  >mpure  the  Ti,:n  i :  -i  squar-  :  :>lot  >  •:  the  m  =  1, 

■  =  ■  >  matched  filter  t  «>’  i  gure  2  .  p  ’  with  th.-  maga  i  t  u.ie  squared  t  he  G  MF 

•>  Figure  2.  Mi.  I’learlv  the  GMF  it  i  !  i  zes  high  spat  in!  t  fqiicn:  i.-s 
pre  ferent  i  i !  !  v.  iiv  complex  con  jnc.it  ion  foil  <wed  hv  inverse  Fourier 
transformation  opera t  ing  on  the  matched  filter  for  the  m  =  1,  =  0  case ,  we 

recreate  tin1  original  scene  (Figure  2.1).  Applying  the  same  sequence  of 
operations  to  the  IMF  leads  to  a  new  scene  which  is  largely  real  but  has 
significant  negative  portions.  The  squared  magnitude  of  that  scene  is  shown 
in  Figure  2.11 .  Here  we  see  even  more  detail  of  those  features  which  are  most 
useful  in  recognizing  the  nine  input  images  and  distinguishing  them  from  white 
no i se . 

Vet  another  wav  of  illustrating  GMF  behavior  is  to  compare  the  OIF 
correlation  plane  with  the  correlation  plane  produced  by  the  m  =  1,  c  =  0 
input  case  (Figures  2.b  and  2.101.  in  Figures  2.12  and  2.13  we  show  the  m  = 
1.1,  o  =  7.b'  case.  Here  the  advantage  of  the  (IMF  in  object  location  is  quite 
obvious.  In  every  case  the  (MF  gives  a  narrow  peak  centered  on  the  airplane 
center  of  mass.  The  matched  filter  always  gives  a  much  broader  peak  and  that 
peak  mav  not  be  centered  properly.  Indeed  we  mav  have  multiple, 
improper  1 v- located ,  broad  peaks. 

Having  exercised  the  GMF  to  achieve  a  broad  case  (all  nine  m,  0 
combinations),  we  then  exercised  it  to  achieve  high  separability.  For  this 
case,  we  kept  ft  =  0  and  used  magnifications  0.8,  1.0,  1.1,  and  1.2.  We 
produced  five  GMF's  each  designed  to  recognize  one  magnification  and 
discriminate  against  the  others.  Here,  matched  filters  do  very  poorly,  see 
Figures  2.3  and  2.4  for  example,  but,  alternative  filters  have  been  suggested 
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Figure  2. 


The  Squared  Magnitude  ut  the  ('nmplex  Image  t.'  Whieh 
the  C.MF  is  "Matched" 
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Figure  2.12.  Ramtek  Display  of  the  Correlation  Plane  Created  Rv  the  Kilt 
Matched  to  e  in  Figure  2.2  When  the  Input  is  C 


bv  Caulfield  and  Malonev.  1  Thus,  we  will  compare  the  (IMF's  with  the 
uaul f ield-Malonev  filters  (CMF's).  The  Cau l f i e 1 d-Ma 1 onev  filters  ( CMF ' s )  are 
linear  combinations  of  filters  matched  to  the  five  differently  magnified 
images.  The  weights  are  so  chosen  that  the  response  to  the  CMF  for  one 
image,  say  the  m  =  1,  image  is  high  but  the  response  to  any  other  m  is  zero 
(at  the  centroid  of  the  image).  In  this  sense  (contrast  between  "true"  and 
"false"),  the  CMF  is  optimum.  figure  2-14  shows  the  correlation  plane 
produced  bv  the  m  =  1  CMF  with  the  m  =  1  input.  Figure  2-15  shows  the 
correlation  plane  with  the  same  CMF  but  m  =  0.8  input.  The  CMF,  though,  is 
defined  only  for  a  finite  set  of  inputs  whereas  the  GMF  is  defineable  for 
infinite  stochastic  sets.  Thus  the  GMF  is  more  versatile.  In  addition,  only 
the  GMF  simultaneously  optimizes  between-c 1  ass  and  within-class  separations. 

2.3  Discussion 

Generalized  matched  filters  appear  to  be  more  flexible  than  matched 
filters  (no  flexibility)  or  Caul f ie Id-Mahoney  filters  (useful  only  for 
increasing  between-c lass  discrimination).  Not  inluded  in  our  study  was  the 
within-class  filters  due  to  Hester  and  Casasent . 4  The  latter  authors  allude 
to  the  possibilitv  of  combining  their  work  with  what  amounts  to  the 
Caul field-Maloney  work  to  achieve  full  flexibility.  Leger  et  al .  ’  also  used 
statistical  methods  to  develop  quite  different  spatial  filters  with  much 
flexibility.  Likewise,  Caulfield  et  al.*’  showed  how  to  weight  matched  filters 
to  have  full  flexibility,  i.e.,  with-class  and  between-c 1  ass  control. 

This  paper  contains  the  first  demonstration  that  GMF's  can  be  derived  for 
very-general  two-dimensional  scenes.  It  shows,  as  well,  the  qui te-general 
property  that  GMF  correlation  peaks  are  much  narrower  than  the  corresponding 
classical  cross  correlation  peaks  produced  bv  matched  filters.  This  results 

from  the  different  senses  in  which  matched  filters  and  GMF's  are 

•  •  ? 

opt imized . 
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Figure  2 .  I  "> .  The  (’.or  re  1  tl  I  •  •  i  Plane  for  the  em  -  1  CMK  with 
m  =  0.8  Input 
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Of  all  the  matched  filter  replacements,  only  the  GMF  absolutely  requires 
computer  recording.  On  the  other  hand,  all  of  the  methods  require  computer 
calculat ion . 

Computer  recording  to  the  required  accuracies  is  beyond  the  capabilites 
of  available  hologram  writers.  We  are,  however,  developing  a  hologram  writer 
which  is  capable  of  writing  10®  and  106  points  to  an  absolute  spatial  accuracy 
of  one  part  in  10®. 7  A  more  subtle  problem  is  the  fact  that  optical  systems 
do  not  necessarily  perform  good  Fourier  transforms.  Indeed  because  MTF ' s 
degrade  off-axis,  real  systems  are  space-var iant .  The  GMF  has  the 
potentiality  of  creating  the  stat i st ical ly-opt imum  space-invariant  filter  for 
that  space-variant  system.  In  addition,  matching  the  spacing  of  the  filter  to 
the  exact  size  spacing  and  orientation  required  for  the  optical  equivalent  of 
our  digital  matching  can  be  a  problem.  Because  the  GMF  emphasises  the  high 
spatial  frequency  domain,  we  might  expect  it  to  be  more  sensitive  to  angular 
alignment  than  the  matched  filter.  This  will  be  the  case  in  many 
circumstances .  In  our  case,  however,  the  exact  opposite  is  true.  The  angular 
insensitivity  to  alignment  of  the  filter  comes  about  from  the  fact  that  the 
GMF  was  designed  to  be  relatively  insensitive  to  rotation  between  the  input 
and  the  filter.  It  is  clear  that  we  could  as  well  design  the  filter  to  be 
insensitive  to  translation  error.  Each  new  insensitivity  we  add  must  be 
expected  to  cost  us  something  in  terms  of  between-c lass  discrimination,  so  we 
have  avoided  adding  that  translation  insensitivity  to  this  particular  filter. 
We  are  currently  at  work  on  the  synthesis  of  computer  generated  GMF's  and  hope 


to 

report 

on  that  in 

the 

near  future. 
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3.  HOLOGRAM  CONSTRUCTION 


3 . 1  Characterizing  the  Need 

Here  we  seek  to  define  the  necessary  characteristics  of  a  hologram  writer 
for  practical  pattern  recognition.  The  primary  quantities  are 

o  N  (the  number  of  pixels  that  can  be  written  in  each  direction), 

o  p  (the  size  of  the  individual  pixel),  and 

o  6p  (the  maximum  absolute  accuracy  of  the  separation  of  pixels 
and  N ) . 

Note  that  perfect  optical  magnification  or  demagni f icat ion  leaves  N 
unchanged  while  changing  p  to  mp  and  Sp  to  m6p.  Accordingly  the  most 
fundamental  properties  of  the  writer  are  N  and  p/6p.  Nevertheless  since 
demagnification  may  lead  to  distortion  which  must  be  precompensated  for  in 
hologram  design,  it  is  important  to  have  5p  small  enough  to  produce  a  usable 
hologram  directly.  It  is  not  necessary  that  6p  be  small  enough  to  introduce  a 
big  splitting  angle,  because  we  can  use  the  object  wavefront  from  the  computer 
generated  hologram  with  a  separately-derived  reference  beam  to  form  an  optical 
hologram  with  any  angle. 

Thus  we  can  arrive  quickly  at  a  minimum  value  for  N.  We  know  that  an 
off-axis  hologram  of  wavefront  M  pixels  across  requires  at  least  4M  pixels. 
Very  good  lenses  give  as  many  as  10,000  pixels,  so  an  off-axis  hologram  for 
use  with  a  good  lens  might  require  N  >  40,000.  Most  lenses  we  consider  for 
optical  pattern  recognition  are  not  that  good,  so  we  do  not  need  such  large 
N.  Nevertheless,  N  =  40,000  is  safe. 

It  is  also  easy  to  set  the  fringe  width  p  with  some  accuracy.  For  the 
greatest  ray  deviation  to  be  8,  we  need  to  be  able  to  write  fringes  of  width 
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X _ 

4  sin  (9/2) 


(13) 


For  small  0  , 


P  “ 


X 

29 


(14) 


For  a  reasonable  splitting  angle  like  0  =  0.1  ("5.7°)  and  X  »  0.6  pm  (red), 
we  have  p  *  3  pm.  Thus  p  ■  3  pm  is  a  reasonable  goal. 

To  represent  an  ideal  wavefront  to  X/Q  accuracy,  we  require  26w  (which 
corresponds  to  the  representation  of  a  full  wave)  to  be  represented  to  1:0 
accuracy,  i.e. 


>  0/2  (15) 

For  a  very  good  hologram  X/20  accuracy  seems  appropriate. 

Therefore  we  start  with  goals  of 

N  ’40,000, 
p  ”3  pm,  and 
<$p  ”0.3  pm. 

Anticipating  what  follows,  the  hologram  writer  we  produced  has 

N  -  14,000, 
p  ■  5  pm,  and 
6p  -  0.15  pm. 

Note  that  it  is  within  a  factor  of  2  to  3  of  ideality  on  all  counts  and  better 
than  need  be  on  6p.  As  the  goals  themselves  are  inexact,  we  conclude  that  the 


3-2 


hologram  writer  we  have  constructed  will  be  adequate  for  all  or  almost  all 
future  needs  for  computer  generated  pattern  recognition  holograms. 

3 . 2  Consequences  of  the  Hologram  Writer  for  Hologram  Encoding 

The  capabilities  of  the  hologram  writer  change  the  way  we  think  of 
hologram  encoding.  Because  Np/6p  =  (1.6  x  10+4  x  5)/0.15  =  5  x  105,  as 
opposed  to  the  more  common  10 3,  the  "phase  quantization"  problem  which  has 
plagued  prior  computer  holograms  is  effectively  negligible.  Because  N  is 
substantially  larger  than  before,  many  more  shades  of  gray  (proportional  to 
N  )  can  be  recorded  for  a  given  number  of  object  wavefront  pixels  than  ever 
before.  Consider  one  of  the  prior-art  CRT  writers  with  p/6p  ..1  at 
smallest  p  and  N  =  2000.  If  we  want  to  represent  400  object  wavefront  pixels, 
we  have  a  4  x  4  cell  for  each.  The  cell  center  can  have  one  of  four  possible 
positions  yielding  A/4  accuracy.  The  cell  power  transmission  can  vary  from  0 
to  1  in  16  steps.  With  our  writer  we  have  a  32  x  32  cell  for  each  or  1024 
shades  of  gray.  The  phase  accuracy  is  X/(32  x  p/ 6p)  »  A/100.  In  the  prior 
art  A/1000  accuracy  would  require  an  M  *  2  hologram  (a  2  pixel  x  2  pixel 
object  wavefront).  Likewise  1024  shades  of  gray  (32  x  32  cells)  would  allow 
only  a  64  pixel  by  64  pixel  object  wavefront.  Clearly  the  new  hologram  writer 
obviates  many  of  the  tradeoff  needs  necessitated  by  previous  writers. 

Now  that  computer  holography  is  a  serious  alternative,  we  must  study  some 
practical  issues  which  were  previously  ignored.  An  example  of  such  a 
consideration  is  lens  imperfection.  Perhaps  the  easiest  way  to  map  this  is  to 
use  the  actual  optical  set  up  to  Fourier  analyze  a  square  input.  Locations  of 
the  peaks  and  valleys  of  the  power  spectrum  (uniformly  -  spaced  in  the  ideal 
case)  will  not  be  uniform.  This  permits  us  a  peculiar  sort  of  "distortion" 
measurement.  We  then  determine  whether  or  not  this  "distortion"  is  invariant 
to: 
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o 


Input  position  (space  invariance)  and 


o  Size  of  the  square  (should  we  call  this  "distortion  invariant 
distortion"?) . 

We  can  then  check  to  see  if  the  power  readings  follow  the  ideal  "sine  squared" 
pattern  and  whether  the  power  pattern  is  space  invariant  or  "distortion 
invariant."  Armed  with  these  nonideal  properties  we  can  design  the  space 
invariant  hologram  optimum  in  some  stochastic  sense  (Bayes,  miniraax,  etc.)  for 
the  system. 


3.3  Basic  Scheme  of  the  Holowriter 


The  basic  scheme  is  to: 

1.  Use  an  ultra-low  distortion  CRT  (greatly  demagnified)  to  write  a 
small,  square  portion  (cell)  of  the  hologram  onto  a  photographic 
emulsion  on  an  ultraflat  photographic  plate, 

2.  Translate  the  plate  to  write  the  next  cell, 

3.  Measure  the  translation  to  an  absolute  accuracy  of  0.1  pm 
inter ferometr ically , 

4.  Write  the  next  cell  on  the  CRT  face  with  position  chosen  to 
compensate  for  the  unpredictable  (to  within  a  few  microns)  stopping 
position  of  the  stage. 


Of  course  such  a  system  must  be  computer  controlled.  Achieving 
meaningful  0.2  im  accuracies  requires  extremely  careful  engineering  as  well  as 
the  use  of  a  temperature  ontrolled  room. 

Figure  3.1  shows  the  overall  system  layout  schematically.  At  least  a 
dozen  factors  converge  to  make  improvements  beyond  this  performance  level  very 
difficult  and  very  expensive.  A  few  of  those  factors  follow: 

i.  The  special  right-angle  mirror  required  for  the  2-D,  orthogonal 
interferometry  can  not  be  made  much  longer  without  itself 
introducing  0.1  um  errors, 
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2.  The  writing  time  for  a  10^  x  lO*’  hologram  is  approaching  a  full 
working  day  and  system  stability  for  longer  exposures  is  doubtful 
(to  say  nothing  of  latent  image  failure), 

3.  Only  one  commercial  photographic  plate  exceeds  minimum  requirements 
on  sensitivity,  resolution,  and  flatness  even  with  the  present 
design , 

4.  The  demagnification  is  already  so  great  that  it,  not  depth  of  focus, 
drives  us  to  use  state-of-the-art  controllers  of  lensplate  distance, 

5.  Temperature  gradients  (beyond  those  already  sensed  and  compensated 
for)  as  well  as  temperature  drift  limit  mechanical  measurement 
accuracies  to  no  better  than  about  0.15  \xa, 

6.  Temperature  variations  within  the  room  limit  the  CRT  demagnification 
to  roughly  its  present  value. 
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4.  DESIGN  CONSIDERATIONS 


4 . 1  Introduction 

Designing  the  holowriter  required  the  careful  selection  of  components  and 
equallv-careful  attention  to  the  environment  in  which  those  components  must 
function.  In  this  section  we  describe  those  considerations. 


4 . 2  Temperature  of  Room 

The  primary  effect  of  changing  room  temperature  is  to  change  the  physical 
dimensions  of  all  of  the  metal  components,  bases,  and  supports.  Listed  below 
are  some  of  the  ways  we  can  expect  those  dimension  changes  to  cause  problems: 

o  The  demagnification  factor  may  vary  because  the  object  distance  and 
image  distance  (relative  to  the  principal  planes  of  the  microscope) 
may  change  as  the  basic  structure  shrinks  or  expands, 

o  The  lateral  position  of  the  optical  axis  of  the  microscope  may  change 
relative  to  the  CRT  and/or  relative  to  the  stage,  and 

o  The  separation  between  the  interferometer  and  the  mirror  it  is 
measuring  may  be  due  partially  to  mirror  motion  and  partiallv  to 
temperature-induced  changes  in  system  (e.g.  baseplate  dimensions). 

There  are  only  four  ways  to  combat  this  problem:  design  in  temperature 
compensation,  minimize  temperature  change,  use  non-expanding  materials,  and 
measure  temperature  and  correct  for  it. 

Compensation  "constancy"  was  achieved  by  building  the  system  in  a 
temperature  controlled  room.  We  also  kept  most  heat  sources  (except  the  CRT) 
outside  that  room.  There  the  questions  are:  how  stable  is  the  temperature 
over  the  anticipated  writing  times  and  how  big  are  the  temperature  gradients? 
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Low  expansion  materials  (e.g.  special  ceramics)  were  seriously  considered 
for  the  baseplate  and  other  major  structural  components.  Cost  and  convenience 
ruled  this  out. 

Measure-and-compensate  mechods  are  used  to  correct  the  interferometer 
readings  for  the  effects  of  temperature  gradients. 

When  we  combine  the  effects  of  all  of  the  trade-offs  and  compensat ions  we 
come  Co  the  conclusion  ChaC  the  later  accuracy,  <$x,  depends  on  the  total 
temperature  excursion  At  during  the  total  writing  time,  tW)  according  to 

6x  *  a  ATn(cw)  .  (16) 

Here  the  proportionality  constant  a  is  about  0.1  pm/'C.  Time  enters  in  as  a 
further  drift  problem.  Thus  there  can  be  a  drift  in  the  reference  to  which 
the  room  temperature  is  compared.  We  believe  that  this  is  not  significant 
during  any  one  working  day. 

In  a  future  and  more-expensive  system  we  could  obtain  significant 
improvements  by: 

o  A  temperature  controller  with  greater  stability  and 

o  Use  of  low  expansion  ceramics  wherever  possible.  This  could  lead  to 
roughly  a  factor  of  two  improvement. 

4 . 3  Stage  Motion 

As  the  stage  micrometers  are  driven,  the  stage  itself  undergoes  periodic 
up  and  down  motion  from  the  gear  drive.  This  motion  is  certified  by  the 
manufacturer  to  be  1  I  pm.  The  primary  effect  is  not  on  focus  (our  depth  of 
focus  is  around  25  pm)  but  on  the  deraagnif icat ion  ratio 

m  -  Sq/Si  ,  (17) 
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where  S0  =  object  distance  and  Si  3  image  distance  for  the  demagnif ying 
microscope.  The  relative  error  in  m  is 


5m 

m 


where  5S0  and  ^Si  are  the  errors  in  So  and  Si.  Stage  motion  of 
<5Si  =  I  pm  =  lO-4  cm  leads  to  a  relative  magnification  error 

-3 

Sra/m  *  5Sj/Sj  *  5  x  10 


(18) 


(19) 


This  gives  Sm/m  =  5  x  10“3  for  Sj  =  1.9  x  10-2  cm.  Note  that  changing  m 
leads  to  a  change  in  the  fringe  position  because  we  are  writing  a  field  of 
nominal  site  200  v*n.  Thus  6m/m  =  5  x  10“ 3  leads  to  a  maximum  error  of  1  pm  in 
fringe  resolution. 

When  we  attempt  to  drive  the  stage  to  a  new  position,  we  arrive  at  a 
location  xSjy„  which  differs  from  the  target  position  xt,yt  by 
unpredictable  errors  6x  and  6y.  That  is 


x  =  x  +  6x  (20) 

s  t 


and 


yt  +  fiy 


(21  ) 


These  errors  are  the  settling  errors  certified  to  be  roughly  10  pm.  Our 
approach  to  dealing  with  these  errors  is  to 


o  Measure  xs ,  Xt, 
o  Calculate  6x,  5y,  and 
o  Write  the  CRT  image 
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laterally  displaced  from  the  center  by  5x/m,  6y/m.  The  demagnified  image  is 

then  automatically  centered  on  x  ,y  . 

t  t 

4.4  Photographic  Plate  Choice 

The  photographic  plate  roust  be  extremely  flat  to  achieve  good  resolution 
because  of  the  previously-noted  demagnification  effect.  Kodak  offers  the 
flatest  available  photographic  plates  (Microf lat®) .  These  have  a  flatness  of 
6S  =  4  tin  over  the  full  usable  aperture  of  7  cm  and  thus  limit  writing 
accuracy  on  fringe  width  due  to  the  plate  (as  explained  in  Section  3.3) 
to  4  van. 

In  photographic  emulsions  big  grains  are  used  to  gain  high  speed. 

Because  we  want  high  resolution  (-5  pm)  we  must  expect  low  speed.  Lower 
speed  means  longer  exposures.  Longer  exposures  means  less  resolution. 

Clearly  there  is  a  narrow  "window"  in  the  set  of  all  emulsions  which  satisfies 
a  mutual ly-consistent ,  useful  tradeoff  between  speed  and  resolution.  Of  the 
few  potentially-interesting  emulsions,  only  one  was  available  on  Microflat 
plate.  We  conclude  that  this  was  the  only  commercially-available  photographic 
plate  suitable  for  use  with  this  hologram  writer. 

4.5  CRT  Choice 

The  CRT  choice  is  dominated  by  considerations  of  speed  and  distortion. 

To  avoid  writing  the  hologram  a  point  at  a  time  we  chose  to  write  it  a  square 
cell  at  a  time.  The  speed  improvement  factor  in  so  doing  is  the  number  of 
points  in  the  cell.  The  ultimate  speed,  of  course,  would  come  from  writing 
the  whole  hologram  at  once.  For  the  number  of  points  and  for  the  positional 
accuracy  we  require,  this  ultimate  CRT  does  not  exist.  The  question  then  is 
how  large  a  cell  (in  terms  of  number  of  resolvable  elements)  can  we  write  to 
the  required  accuracy.  Our  original  goal  was  to  write  fringes  to  A/100 
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accuracy.  With  a  fringe  period  of  10  \sa  this  requires  0.1  p  on  the 

hologram.  On  the  CRT  it  requires  p/50  accuracy  across  the  cell,  where  p  is 

the  width  of  a  single  point.  That  is,  we  need  to  confine  the  cell  to  that 

region  of  the  CRT  giving  less  than  2%  distortion.  For  a  band  CRT  the  cell  may 

be  very  small.  For  the  best  CRT,  it  may  be  very  large.  For  the  good  and 

affordable  CRT  we  chose  the  Infodex  PD1200M34  with  a  12  cm  diameter  RCA  C82200 
ESI  tube.  The  usable  area  was  5  cm  x  5  cm  for  our  purposes. 

The  image  must  be  "on"  long  enough  to  expose  the  emulsion.  We  hope  it  is 

bright  enough  so  that  retrace  is  not  required  to  achieve  the  desired  exposure 
level.  The  phosphor  decay  time  would  be  less  than  the  time  needed  to  step 
the  stage  and  let  it  settle  to  with  0.1  in.  For  our  system  that  time  is 
0.17  seconds.  The  decay  time  for  the  RCA  tube  phosphor  is  0.001  seconds. 

The  chief  obstacle  to  better  CRT  performance  (bigger  cells)  is  cost.  The 

next  best  CRT  was  roughly  $5,000  more  than  the  one  we  bought.  In  a  program  in 

which  cost  was  a  prime  consideration,  this  was  not  a  serious  option. 

4.6  Mirror  Design 

The  right  angle  mirror  attached  to  the  stage  to  allow  stage  x-y  position 
monitoring  must  be  locally-flat  to  allow  the  interferometer  to  function.  In 
addition  it  must  be  accurately  90°  in  order  to  not  couple  x  and  v 
measurements.  Suppose  it  is  accurately  90°,  we  can  align  the  mirror  at 
x  =  x0  and  adjust  the  mirror  alignment  so  that  the  x  indication  is  constant 
to  better  than  0.1  pm  regardless  of  y.  The  test  is  to  move  to  a  much 
different  value  of  x  and  scan  y  there.  It  should  again  lead  to  "constant"  x. 
If  the  angle  is  not  90°  but  90°  +  69,  then  the  maximum  "coupling  error"  is 

ce  =  l  68  (22) 
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where  l  is  the  maximum  scan  distance  (around  8  cm  in  this  system).  Thus  we 
would  like  to  have  ce  <  0.1  yn  or 


66  <  10"5  cm/8  cm  =  1.25  x  10~6  .  (23) 

Microradian  accuracy  over  8  cm  is  extremely  hard  to  achieve  (or  even  verify). 
Testing  in  the  visible  across  the  whole  face  gives  a  diffraction  limit  of 

A9  -  \/l  «=  0.5  x  10"4/8  "566  .  (24) 

It  follows  that  we  can  not  expect  to  avoid  correcting  x  according  to  y,  i.e. 
calibrating  the  mirror  and  storing  a  simple  correction  factor  for  use  by  the 
computer . 

4 . 7  Vibration 

In  a  system  designed  to  approach  0.1  im  positional  accuracy,  normal 
vibrations  would  be  catastrophic.  We  have  taken  four  precautions  to  minimize 
vibration.  First,  we  have  used  an  air-suspended  table  to  isolate  the  system 
from  normal  slow  vibrations  (people  walking,  machinery  operating,  earth 
settling,  etc.).  Second,  we  have  used  a  granite  table  to  damp  fast 
vibrations.  Third,  we  have  used  only  heavy,  solid  components.  Fourth,  we 
write  holograms  only  after  disturbances  from  stage  motion  and  disturbances 
from  human  operators  have  had  time  to  dampen  out. 

4 . 8  Air  Currents 

The  temperature  controlled  room  continually  circulates  air.  This  should 
maintain  the  desired  temporal  and  spatial  uniformity,  but  it  can  also  cause 


sufficient  turbulence  to  cause  misdirection  of  the  writing  beam.  We  hav<> 
built  turbulence-spoiling  shields  to  remove  this  problem. 

4.9  Interferometer 

The  interferometer  was  one  of  the  easiest  choices.  What  we  sought  was 
c  0.1  (jn  resolution  or  better 

o  reading  speed  c.fficient  to  verify  damping  of  stage  motion  in  a  small 
time . 

The  standard  of  the  industry  is  the  Hewlett-Packard  interferometer  we  chose. 

It  is  more  than  adequate  on  both  counts. 

4 . 10  Radio-Frequency  Interference 

We  anticipated  that  radio  frequency  interference  (RFI)  would  present  a 
problem  for  our  delicate  display  circuitry  and  tried  to  shield  the 
environmentally-controlled  room  against  it.  It  turned  out  that  RFI  also 
affected  the  controlling  PDP  11/23  computer.  After  much  effort,  this  and  a 
line  surge  problem  caused  by  the  room  controls  was  brought  under  sat  is fac tor y 
control . 

4.11  Information  Flow 

We  needed  to  accomplish  two  primary  tasks  by  computer: 

o  Calculation  of  the  hologram  (a  nontrivial  task  in  that  the  writing 
beam  must  be  fed  information  to  0.1  un  accuracy  over  an  8  cm  diameter 
hologram  for  a  total  of  about  5  x  10"  bits),  and 

o  Control  the  writing  (move  the  stage,  sense  stage  position,  write  the 
cell,  decide  next  cell  position,  obtain  cell  data,  etc.). 

In  principle,  one  computer  could  do  both  tasks.  In  practice,  this  proved  far 
too  expensive. 
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To  operate  the  holowriter  at  its  design  speed  while  making  fringe 
calculations  on-the-fly,  Aerodyne's  PRIME  computer  would  have  to  be  devoted 
full  time  to  the  task  (all  other  users  banned  for  th?  entire  writing  time). 
The  financial  burden  to  the  program  of  this  operation  would  have  been 
immense.  We  chose,  instead,  to  split  the  tasks  between  two  computers  by 
providing  a  dedicated  PDP  11/23  to  control  the  holowriter.  The  primary  cost 
(in  time,  inventiveness,  and  dollars)  this  approach  caused  was  that  an 
interface  between  the  computers  had  to  be  designed  and  built  and  a  100  foot 
special  cable  had  to  be  designed  and  installed  to  connect  the  computers.  A 
simplified  task  division  is  shown  in  Figure  4.1. 


II 


5.  SYSTEM  DESCRIPTION 


5 . 1  I n t rod  uc t i on 

Having  described  the  system  philosophy,  the  system  components,  and  the 
system  design  parameters,  we  turn  now  to  the  description  of  the  system  as  a 
whole.  In  a  sense,  the  system  was  designed  once  the  components  were 
selected.  Unfortunately  most  of  the  hard  work  went  into  assembly  and 
integrat ion . 

5 . 2  System  Appearance 

Figure  5 . 1  shows  the  CRT,  mount  optics,  stepping  motors,  and 
interferometer  on  the  granite,  air-suspended  table  in  the  envi ronment a  1 1 y- 
eontrolled  room.  The  room  environment  is  set,  controlled,  and  recorded  with 
the  equipment  shown  in  Figure  5.2.  Two  entries  to  the  room  are  locked  during 
hologram  writing  (Figure  5.3).  This  room  had  to  have  all  light  leaks  plugged 
and  all  control  lights  covered  because  the  plate  is  panchromatic  and  is 
exposed  for  a  multihour  period.  All  manual  man i pu 1  a t i ons  after  the  plate  is 
removed  from  the  box  are  done  in  darkness  using  an  infrared  viewer.  Outside 
the  second  door  of  the  primary  double-chambered  entry  wav  (Figure  5.4)  is  the 
PDP  11/23  and  associated  I/O  equipment.  The  operator  controls  all  operations 
from  here.  Down  the  hall  is  the  PRIME  computer  (Figure  5.5)  used  for  the 
fringe  calculations  and  the  interface  based  on  the  PRIME/RAMTEC  interface 
(Figure  5.6)  which  was  already  designed  for  rapid  data  transfer  from  the 
PRIME. 

5 . 3  Alignment,  Magnification,  and  Focus 

Plates  are  inserted  in  the  dark  (as  noted  above)  and  kinematically 
located  within  the  plate  holder.  The  location  of  the  first  fringe  is  selected 
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igure  t.<'.  The  K.imt  ok  Displav  (a)  is  interfaced  tn  the  Prime  Through 
Special  Boards  Designed  and  Built  at  ARI.  A  single  panel 
lit  additional  electronics  (hi,  converts  that  interface  to 
the  required  high  data  rate  Prime  I’DP  ll/HB  interface. 


by  the  starting  point  on  the  micrometers  on  the  stage  motion  transducers. 

The  magni f icat ion  can  be  "trimmed"  electronically.  By  exposing  square  cells 
which  ought  to  abutt  perfectly  while  changing  the  magnification  electronically 
between  cell  exposures,  we  can  achieve  agreement  between  the  cell  size  on  the 
plate  and  the  interferometrical ly-cont rol led  stage  motion.  Thus  the  scale  of 
the  cell  can  be  accurate  to  interferometric  accuracy.  Focus  can  be  obtained 
through  direct  observation  of  backscatter  through  the  pellicle  beam  splitter. 
It  has  proved  satisfactory  in  practice  to: 

o  Align  the  plate  holder  apparatus  once  optically  so  that  no  focus 
change  occurs  over  the  whole  plate, 

o  Measure  the  air  gauge  gap  indicator  reading,  and 

o  Use  the  air  gauge  to  reset  the 'focus  for  new  plates  (the  microscope 
must  be  elevated  and  reset  in  the  dark  to  insert  a  new  plate). 

5.4  Final  System  Specifications 

Table  5-1  shows  final  system  performance  levels. 

Table  5-1  -  Final  System  Performance  Levels 


Parameter  Value 


Hologram  Size 

7  cm  x  7  cm 

Fringe  Size  (Minimum) 

5  pm 

Absolute  Location  Accuracy 

0.15  pm 

Writing  Time  (Full  Format) 

12.5  hours 

5.5  Time  Budget  For  One  Cell 

Table  5-2  shows  the  time  required  for  the  various  operations. 
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Table  5-2  -  Time  of  the  Writing  of  a  Single  Cell 


Operat i on 

Time  (milliseconds) 

Step 

170 

Expose  (2  traces) 

100 

TOTAL 

270 

5 . 6  Sample  Rulings 


Figure  5.7  shows  a  7  cm  x  7  cm  Ronchi  ruling  with  10 
Figure  5.8  shows  a  magnified  area  covering  several  cells, 
accuracy  and  quality  now  achievable. 


pm  lines. 

This  illustrates  the 
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APPENDIX  A 

DETAILED  PROCEDURE  FOR  CALCULATING  THE  GMF 


We  describe  here  the  detailed  procedure  which  would  allow  any  interested 
reader  to  generate  his  own  GMF's. 

We  chose  to  work  with  a  128  x  128  grid,  but  this  choice  was  arbitrary. 

We  could  go  larger  or  smaller.  Binary  numbers  are  necessary  for  doing  fast 
Fourier  transforms  (FFT's). 

The  GMF  algorithm*  works  with  sampled  values  of  the  FFT  at  all  points  in 
the  filter  plane.  Our  program  allows  us  to  work  with  either  the  actual 
sampled  points  or  with  a  mean  and  Gaussian  variance  calculated  from  them.  We 
chose  the  latter  approach.  Thus,  for  each  of  128  x  128  points,  we  calculate 
and  store  the  mean  and  variance  of  FFT’s  at  that  point  the  associated  various 
inputs  scenes  (in  this  case,  rotated  and  magnified  aircraft). 

The  FFT's  will  have  both  real  and  imaginary  parts.  We  treat  them  as 
independent  128  x  128  matrices.  The  real  matrix  generates  the  real  part  of  the 
GMF  and  the  imaginary  matrix  generates  the  imaginary  part. 

We  do  not  need  to  calculate  the  2  x  128  x  128  GMF  values  simultaneously. 
We  can  break  them  up  arbitrarily  and  get  the  same  result  because  we  are 
assuming  (precisely  for  this  immense  calculational  advantage)  that  values  at 
each  of  these  points  are  independent  of  the  values  at  any  other.  We  chose  to 
run  the  256  rows  (of  128  points  each)  sequentially  and  store  the  resulting 
values  of  the  "weights"  s  coefficients  in  the  GMF. 

Actually,  what  we  calculate  in  this  way  is  not  the  GMF  but  its  complex 
conjugate.  Thus,  to  obtain  the  correlation  pattern  for  any  input  scene,  we: 


*The  Actual  algorithm  is  Program  360A-CM-03X  from  the  IBM  SYSTEM/360 
Scientific  Subroutine  Package  (1970),  p.  52. 


(1)  Perform  an  FFT  (128  x  128)  on  that  scene, 

(2)  Multiply  the  FFT  by  the  GMF,  and 

(3)  Inverse  FFT  the  product. 

Thus,  the  OIF  Is  used  exactly  as  the  matched  filter  is. 

If  we  write 

GMF  +  FFT*  [p(x,Y>] 

then 

p(x,y )  =•  FFT_1  [Off*] 

Here,  p(x,y)  is  the  pattern  to  which  the  Off  is  matched.  As  shown  in 
Figure  2-11,  |p(x,y)|  looks  somewhat  like  an  edge-enhanced,  blurred  image  of 
an  airplane. 
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APPENDIX  B 


Thermal  Expansion  Compensation 

The  attached  analysis  is  a  summary  of  the  final  design  for  the  subject 
equipment . 

The  materials  and  geometry  were  chosen  for  thermal  expansion  symmetry  at 
the  center  of  the  carriage  stroke.  Plus  and  minus  errors  will  occur  at  the 
extremes  of  carriage  stroke.  Also,  the  compensation  was  based  on  uniform  bulk 
temperature  changes  of  the  structure. 

For  the  foregoing  reasons,  two  of  the  the  Hewlett-Packard  ±  0.05°F 
sensors  (#10563AOPT-001 )  were  located  to  monitor  the  temperature  of  the  base 
and  the  aluminum  spacer. 
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Motion  of  elements  is 
referenced  from  zero  datum 
at  optical  center  line. 

6q  *  object  displacement  (CRT) 
6^  =  lens  displacement 

6  *  interferometer  displace- 

x  ment 


The  displacement  error  (5^)  of  the  image  is: 


For  the  CRT 


6i  (1  +  $ 


6^  -  1.007  6^  For  the  lens 


6.-6 
i  x 


For  the  interferometer 


Thus,  thermal  motion  of  the  interferometer  and  lens  is  critical. 
The  CRT  is  no  problem. 


Lens  Motion  (6^) 


For  uniform,  bulk  t  of  the  structure,  we  can  select  the  lens  spacer  to 

compensate  for  6  -  6  .  From  diagram: 

x  y 

6,  “6.  *  6,  sin  45°  and  we  want 

Lx  Ly  L 

&  **  6,  for  compensation  so 

x  Lx 

6  35  4 

6  *=  ■  =  ■  *■  -  “  50  microinches/  F  required,  and 

L  • /U /  • /w/ 

<$L  *  At(L2a2  +  L3a3  -  2.^)  = 

a2  = 
a3  = 

6l  *  At  [2.75  x  10.6  +  13  *3  -  (*3  +  2.75)]  22  * 

6l  =  At  [29.2  +  13&3  -  5.5£3  -  15]  ll  “ 

fiL  =  At  [14  +  7.5£3]  and  we  want  <5^  ■  50  microinches/°F,  so 

*3  *  507~5'1'4  *  4*8"  and  £1  =  7,55  inehes 


5.5  x  10_O/°F 
10.6  x  10~6  for  lens 

13.10~6/°F  for  alum 

2.75  fixed 

*2  +  S  "  *3  +  2,75 


APPENDIX  C 


Photographic  Plate 

The  high  speed  holographic  plates  from  Kodak  (Type  131-02)  on  Microflat 
plates  were  the  only  satisfactory  plates.  Relevant  catalog  data  is  attached. 
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A  noiogram  .3  termed  by  recording  the  nterference  pattern  of 
conerem  gnt  waves.  as  either  amplitude  or  pnase  vanat  ons  in  a 
recording  medium  Typ.caiiy.  one  of  tne  waves  s  either  a  plane 
wave  or  a  divergent  spnencai  wave  ana  it  is  Known  as  me 
reference  peam  The  other  wave  usually  reflecteo  from  an  object 
jnoer  study  s  reterreo  to  as  me  signal  oeam  When  tne  pro¬ 
cesses  hoiogratn  ,s  reuiummatea  Oy  its  original  reference  oeam. 
its  original  signal  oeam  will  oe  reconstructed  m  ooth  amplitude 
ana  pnase  Decause  of  a  hologram  s  unique  reversal  property 
Conversely,  when  'eillummated  by  tne  more  complicated  signal 
oeam.  the  hologram  wai  reconstruct  a  somewnat  modified  but 
similar  reference  beam  This  is  used  m  complex  spatial  filtering 

ideally  holographic  materials  snoud  have  the  following  proper' 


t,es  resolution  greater  than  tne  spatia.  frequency  me 
noiogram  s  interference  pattern,  very  low  grain  noise,  dimen¬ 
sional  staoiiny.  realtime  development  wnnout  processing. 
erasaDiiity  tor  reuse,  no  :n, cur, ess  limitations,  and  mgr,  sensitivity 
at  popular  wavelengths  A  number  ot  materials  Cbmbine  some 
these  properties,  but  none  nas  them  ail 

When  selecting  the  material  most  suitable  to  your  needs  me  two 
most  important  considerations  are  resolution  and  sensitive, 
However,  n  particular  applications  the  support  mg  substrate 
emulsion  thickness  granularity,  maximum  contrast  density,  and 
cost  may  al:  be  important  The  chans  on  tne  opposing  page  w,. 
help  you  select  the  proper  material  tor  your  application  Turn 
prices  are  on  page  94 


SUBSTRATES:  The  primary  suostrates  used  to  Support 
noiograpmc  emulsions  are  glass  plates,  acetate  sheeting,  and 
polyester  sneetmg  Glass  piates.  being  tn.y  most  rigid  and  stable 
can  be  neid  m  a  variety  or  holders  and  oe  precisely  repositioned 
attar  development  Holographic  films  in  roll  or  sheet  format  are 
ess  expensive.  Out  difficult  to  noio  stationary  tor  long  exposures 
Polyester  oases  ,sucn  as  the  Kodak  ESTAR)  are  strong  ana 
dimer.s.onaiiy  stable  out.  unfortunately,  they  exmbit  oiretrigence 
wh.ch  may  cose  problems  in  some  aopi. cations  suen  as  coherent 
optical  data  processing.  For  these  applications,  acetate  base 
•ilms  are  recommended  if  an  emulsion  is  not  available  on  giass 
substrates,  we  advise  tnat  you  use  a  vacuum  film  holder,  or  simply 
sandwich  the  film  between  two  glass  plates  - 

EXPOSURE:  A  primary  parameter  used  to  compare  holographic 
film  s  the  exposure  energy  (ergs, cm2  ;  required  to  produce  a 
ho'ogram  with  a  neutral  density  (01  of  1  0  Generalfy.  in  amplitude 
noiograms.  maximum  diffraction  eHiciency  occurs  between 
0  =0  6  and  0  8.  However,  m  processes  wnere  phase  holograms 
are  generated  by  bleaching  the  image,  a  slightly  nigne:  exposure 
ievei  ,s  required.  Obtaining  the  best  exposure  is  normally  a  trial 
and  error  series  of  tests  However  with  tne  use  of  the  NRC  Model 
EM-1  Exposure  Msek  (pg.  80)  six  exposures  are  obtained 
simultaneously  on  a  single  plate  and  thus  saves  t.me  and  money 
n  selecting  the  optimum,  exposure  tor  a  particular  application, 

RESOLUTION:  A  commonly  used  parameter  to  compare  the 
-^solution  of  pnotographic  film  is  termed  the  'resolving  power". 
Tn  s  value  in  unes/mm)  is  obtained  by  pnotographing  Target  Ob- 
,ects  with  Contrasts  (TOC)  of  either  1.000:1  (high)  or  16  1  (low). 
The  recording  of  holographic  fringes  13,  however,  a  different  pro¬ 
cess  11  has  been  found  that  holograms  were  obtained  even  when 
the  fringes  generated  by  the  signal  and  reference  beams  have 
spacing*  two  to  three  times  smaller  than  that  of  the  photographic 
-esoivmg  power  tor  a  given  film.  Nevertheless,  for  the  saKe  of  pro¬ 
viding  o  common  reference  point,  the  photographic  resolving 
power  has  been  used  here  in  the  chart. 

SHRINKAGE:  Because  it  causes  a  grating  spacing  change,  emul¬ 
sion  snnnxage  has  undesirable  effects  For  example.  shrinkage  ot 
t-ansmission  noiogram  emulsions  causes  small  image  distortions 
'esuitmg  n  spurious  fringes  in  -eal-time  interferometry  Unfor¬ 
tunately.  snnnxage  cannot  De  completely  eliminated,  since  voids 
are  always  created  when  unexposed  silver-halide  grains  are 
romoveo  during  fixing  However  some  success  has  been 
reported  filling  these  voids  with  chemicals  Of  all  the  emulsions 
listed  m  the  film  selection  table.  Kodak  SO-173  appears  to  shrink 
tne  east 


posed  grains  become  giobuies  or  filaments  of  metamc  silver,  that 
introduce  a  certain  amount  ot  unaes.reo  scattering.  Tne  gra¬ 
nularity  values  m  tne  film  selection  table  are  1.000  times  the  stan¬ 
dard  deviation  of  density  about  0=1  when  scanned  by  a  microoen- 
sitometer  having  a  circular  aperture  of  either  48  mm  or  6  mm 
However,  because  of  the  very  iow  granularity  ot  holographic 
emulsions,  reliable  measurements  are  often  hindered  by  .nstru- 
ment  noise 

CONTRAST:  Since  most  nongraphic  materials  have  uniformly 
smaii  grains,  photographic  contrast  is  high  When  low  contrast  ■$ 
oesireo.  you  can  use  either  special  developers  or  simply  lower  the 
developer  concentration,  temperature  anc  development  time  Tne 
contrast  values  (7 )  snown  m  the  fi,m  selection  taoie  identify  me 
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GRANULARITY:  Rnotographic  emulsions  are  fine  silver-halide 
crystals  dispersed  in  a  gelatin  binder  After  processing,  the  ex- 
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r or  noiograpruc  nondestructive  testing  ana  similar  real-time  ap-  light  ana  nave  long  sneit  life .  After  mixing  the  ignt. sensitive 

p.ications  the  NRC  n-Senes  Photopoiymer  provioes  rapiO-ac-  quia  is  spreaa  between  two  giass  piates  to  form  an  approximate. > 

cess  recoramg  with  high  resolution  The  two  important  advan-  10  m  m  thick  layer  The  material  is  neic  in  ptace  Dy  surface  tension 

tages  over  conventional  recording  media  are:  (1)  the  exposed  when  placed  vertically  m  a  plate  noiaer 

noiogram  can  De  vieweo  instantaneously  without  processing,  and 

2)  the  cost  per  snot  can  De  low  for  volume  users.  Exposure*  can  SPECTRAL.  RESPONSE  &  SENSITIVITY:  The  photopoiymer  s 

oe  made  with  argon,  neuum-neon.  or  ruOy  laser  light,  depending  available  in  two  spectral  sensitizations 

cn  tne  spectral  sensitization  chosen  PhOTOPOLYMER  h-i  He-Ne  ano  RuPy  laser  wavelengths 

PHOTOPOIYMER  H-2  Argon  laser  at  488  nm  .me 

Photopoiymer  recording  materials  consist  of  a  solution  of  acrylate  The  energy  required  for  exposure  is  4  to  5  mJ,cm^One-fourtn  of 

monomers,  a  catalyst,  and  a  dye  sensitizer.  Sensitized  dye  that  energy  can  oe  provided  Oy  pre-exposing  or  oiasmg  with  an 

molecules,  when  excited  Dy  laser  iignt,  react  with  the  catalyst  to  incoherent  light  source.  After  exposure,  the  hologram  can  oe  tem- 

produce  free  radicals  that  initiate  polymerization  of  the  acrylate  poraruy  fixed  Dy  incoherent  light  or  permanently  fixed  Dy  heating 

monomers  Therefore,  the  liquid  photosensitive  material  becomes  to  f  50°C  for  to  seconds, 

a  plastic  solid  after  exposure,  image  information  s  stored  via 

modulation  of  the  thicx.ness  and  index  of  refraction  of  the  RESOLUTION  AND  IMAGE  QUALITY:  Nonlinearity  of  the  image 

polymer  and  thus  is  classified  as  a  phase  rather  than  an  approaches  that  of  an  .deal  phase  noiogram.  With  a  resolution  of 

amplitude  recording  material  in  comparison  to  Kodak  649F  Detween  2.000  to  3.000  lines.mm.  image  contrast  of  4Q  to  1  can  be 

plates  the  photopoiymer  is  about  two  orders  of  magnitude  less  Obtained,  in  aoaition,  approximately  20%  of  the  peak  diffraction 

sensitive  but  has  up  to  rO  times  higher  diffraction  efficiency  This  efficiency  is  accessible  within  0.1  sec.  Initially,  the  hologram  may 

allows  a  corresponding  increase  in  the  relative  power  in  the  appear  miiky,  but  after  a  short  period  of  time  will  become  clear, 

reference  beam  to  reduce  the  exposure  time  accordingly. 

PRICE:  Available  m  two  sizes:  25  ml - $60.00 

PREPARATION:  The  photopoiymer  is  sh.pped  in  two  containers  100  ml - $150.00 

Prior  to  mixing,  the  .ndividuai  solutions  are  quite  insensitive  to 


NflC  is  a  dealer  in  Doth  Kodak  and  Agfa-Gevaert  products  for 
nolograpny  and  laser  photography  and  is  stocking  certain 
materials  for  your  convenience  All  materials  are  stored  under  the 
proper  conditions  prior  to  shipping.  Suggested  list  prices  are  as 
of  August.  1977  and  subject  to  revision  without  notice. 

Shipment  from  stock,  quickly  to  anywnere  m  the  continental 
U  S.A 


GLASS  PLATE 


Kodak 

649F  2  x  2  in.  x  0.040  in.,  36  per  box  $52.20 

649F  2  x  2  in.  x  0.040  in.,  Backed,  36  per  box  $57.60 

649F  2  x  2  in.  x  0.060  In.,  Precision  36  per  box  $88.20 

649F  2  x  2  in.  x  0.060  in.,  Preciaion  Sacked,  36  per 

box  $91.80 

649F  4  x  5  in.  x  0.040  in.,  36  per  box  $73.80 

649F  4  x  5  in.  x  0.040  m..  Backed,  36  per  box  $81.00 

649F  8  x  10  in.  x  0.040  in.,  12  per  box  $64.20 

12002  2  x  2  m.  x  0.040  in.,  36  per  box  $52.20 

12002  2  x  2  in.  x  0.060  in.,  Preciaion,  36  per  box  $188.20 

12002  4  x  5  in.  x  0.040  in.,  36  per  box  $73.80 

12001  4  x  5  in.  x  0.040  in..  Backed,  36  per  box  $81.00 

12002  4  x  5  in.  x  0.250  in.,  Microflat,  6  per  box  $59.70 

12002  5  x  7  m.  x  0.250  in.,  Microflat,  6  per  box  $80.00 

12002  8  x  10  in.  x  0.040  in.,  12  per  box  $84.20 

13102  2  x  2  in.  x  0.040  in.,  36  per  box  $52.20 

13102  4  x  5  in.  x  0.040  in.,  36  per  box  $73.80 

13101  4  x  5  in.  x  0.040  in.,  Backed,  36  per  box  $81.00 

13102  8  x  10  in.  x  0.040  in.,  12  per  box  $84.20 

1250  2  2  x  2  in.  x  0.040  in.,  38  per  box  $52.20 

12502  4  x  5  in.  x  0.040  in.,  36  per  box  $73.80 

12501  4  x  5  in.  x  0.040  in.,  8acked,  38  per  box  $81.00 

12502  8  x  10  in.  x  0.040  in.,  12  per  box  $64.20 

12501  8  x  10  in.  x  0.040  in.,  Backed,  12  per  box  $70.80 

HRP-2  4  x  5  in.  x  0.250  in.,  Microflat,  6  per  box  $25.00 

HftP-2  5  x  7  in.  x  0.250  in.,  Microflat,  6  pei  box  $30.65 


Agfa-Gevaert 

8E75NAH  4  x  5  in.,  12  per  box  $14.75 
8E75AH  4x5  in.,  Backed,  12  per  box  $16.20 
8E75NAH  8  x  10  m.,  6  per  box  $19.10 

8E75AH  8  x  10  in..  Backed,  6  per  box  $21.00 

10E56NAH  4x5  in..  12  per  box  $14.75 
10E56AH  4x5  jn.,  Sacked,  12  per  box  $16.20 

10E56AH  8  x  10  m.,  Backed,  6  per  box  $21.00 

10E75NAH  4x5  in.,  12  per  box  $14.75 
10E75AH  4x5  in.,  Backed,  12  per  box  $16.20 

10E75NAH  8  x  10  in.,  6  per  box  $19.10 
10E75AH  8  x  10  in.,  Backed,  6  per  box  $21.00 


SHEET  AND  ROLL  FILM 


Kodak 

649F  35  mm  x  100  ft,  A58  Acetate  Base,  SP414, 

Unperf  $70.80 

SO-1 15  4x5  in.,  50  per  box  $12.30 

SO-173  4  x  5  in.,  25  per  box  $25.20 

SO- 173  35  mm  x  150  ft,  SP417,  Perf  $52.45 

SO-173  9.5  in.  x  200  ft,  SP983,  Unperf  $466.30 

SO-173  12  in.  x  150  ft,  SP983,  Unperf  $442.80 

SO-253  4  x5  in.,  25  per  box  $25.20 

SO-253  35  mm  x  150  ft,  SP41 7,  Perf  $52.42 

SO-253  70  mm  x  150  ft,  SP474,  Unperf  $103.40 

SO-424  4  x  5  in.,  Minlcard  II,  50  per  box  $12.00 

SO-424  36  mm  x  100  ft.  Minicard  II,  SP417,  Perf  $14.00 

2474  4  x5  in.,  Linagraph  Shellburst,  25  per  box  $5.95 

2476  35  mm  x  150  ft,  Linagraph  Shellburst,  Backed, 

SGE417,  Perf  $20.05 

Agfa  Gevaert 

10E75NAH  70  mm  x  100  ft,  Unperf  $44.95 
10E56AH  70  mm  x  100  ft,  Unperf,  $44.95 


APPENDIX  D 


Hologram  Display  Generator 

Attached  is  the  manual  on  the  hologram  display  generator  created 
specially  for  this  program. 
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1.0  INTRODUCTION 


The  Hologram  Display  Generator  is  part  of  a  larger  system  designed  to 

produce  large  format,  high  resolution  holograms  on  glass  photographic 
plates  by  a  step  and  repeat  method.  The  purpose  of  the  HDG  is  to 
generate  single  frame  images  on  a  CRT  which  are  imaged  on  the 
photographic  plate.  Data  describing  the  image  are  transferred  in 
compressed  form  from  a  LSI  11  computer. 

This  manual  includes  a  description  of  the  requirements  imposed  on  the 
HDG  by  system  constraints,  a  functional  description  of  the  HDG 

hardware,  and  an  explanation  of  the  microcode  used  to  satisfy  these 
requirements. 

A  single  hologram  consists  of  about  300  by  300  frames.  Each  frame  has 
a  resolution  of  2000  pixels  by  800  scan  lines  resulting  in  a  total 
information  content  of  1.4E11  pixels.  This  fact  dictates  3  important 
characteristics  of  the  HDG:  high  speed,  efficient  data  encoding,  and 
reliability. 
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2.0  REQUIREMENTS 

The  function  of  the  Hologram  Display  Generator  (HDG)  is  to  accept 
compressed  data  describing  a  given  type  of  hologram  from  the  system 
LSI  11  computer  and  to  generate  the  requisite  signals  to  display  the 

information  in  expanded  form  on  an  Infodex  high  resolution  display 
system. 

The  image  to  be  generated  on  the  display  comprises  a  series  of  bars  of 
varying  width  and  constant  brightness  which  are  parallel  to  one  axis 
of  the  display.  The  image  is  formed  in  a  raster  and  the  bars  are 
defined  in  terms  of  their  widths  in  one  direction  (X)  and  in  number  of 
raster  lines  in  the  orthogonal  direction.  The  image  is  formed  in  a  50 
mm  by  50  mm  square  approximately  in  the  center  of  the  5  inch  diameter 
CRT.  For  purposes  of  this  system,  the  format  represents  2000  pixels 
in  the  X  direction  and  800  scan  lines  in  Y.  The  image  can  be  moved  +/- 
10  percent  in  any  direction,  so  the  useful  (or  used)  portion  of  the 
screen  is  60  mm  by  60  mm  corresponding  to  2400  pixels  by  960  scan 
lines. 

The  signals  required  to  drive  the  display  are  X  and  Y  deflection 
voltages  and  a  binary  unblanking  signal.  The  information  transferred 
from  the  LSI  11  computer  includes  digital  words  describing  the 
position  of  the  origin  of  the  image  (XOFFSET ,  YOFFSET) ,  the  size  of 
the  image  on  the  CRT  (XMAG,  YMAG),  the  number  of  scan  lines  in  each 
group  of  bars  (LRPT ) ,  a  control/status  word  (CSR),  and  a  block  of  data 
of  up  to  1024  words  describing  bar  widths. 

The  time  required  for  a  single  scan  line  is  52  microseconds  which 
includes  30  microseconds  of  active  scan  time,  6  microseconds  for  the 
sweep  to  stabilize,  and  16  microseconds  for  retrace.  The  time  to 
generate  a  single  cell  (one  complete  raster),  if  ail  800  lines  are 
used,  is  41.6  milliseconds. 

The  X  and  Y  deflection  amplifiers  in  the  display  are  identical  and 
have  input  impedances  of  about  1000  ohms  and  deflection  factors  of 
about  .5  inches  per  volt  at  the  input.  The  unblanking  amplifier  has 
an  input  impedance  of  50  ohms  and  a  required  voltage  swing  of  0  to 
-3.5  volts  where  0  volts  corresponds  to  a  blanked  condition.  The 
display  has  a  front  panel  brightness  control,  sc  brightness  need  not 
(indeed,  cannot)  be  controlled  by  the  HDG. 

The  display  is  located  in  a  controlled  environment  chamber.  Because 
of  the  speed  and  precision  of  the  signals  between  them,  the  HDG  is 
located  adjacent  to  the  display.  The  computer  is  located  outside  the 
controlled  chamber  about  20  cable  feet  away. 


\ 
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3.C  DETAILED  DESCRIPTION 


The  major  elements  of  the  HDG  are  shown  in  the  block  diagram  in  Figure 
1.  Interface  elements  between  the  HDG  and  the  Infodex  display  and  the 
LSI  11  computer  are  evident  in  the  block  diagram.  The  memory  is  used 
to  store  the  block  of  data  describing  bar  widths.  These  descriptors 
are  converted  to  an  unblanking  signal  by  the  counter  and  delay 
circuits  and  the  number  of  scan  lines  in  a  bar  is  counted  in  the  line 
repeat  counter.  The  various  HDG  circuits  are  controlled  by  a 
microcontroller  in  accordance  with  a  ROM  stored  program. 


3.1  LSI  11  Interface 

The  interface  between  the  KDG  and  the  LSI  1 1  computer  includes  a 
Programmed  I/O  port  for  transferring  the  single  word  information  and  a 
DMA  port  for  transferring  the  block  of  words  describing  bar  widths. 

The  PIO  interface  is  split  between  a  quad  width  wire  wrap  board 
located  in  the  computer  and  the  main  HDG  circuit  board.  The  two  are 

connected  by  two  20  foot  ribbon  cables.  The  base  address  and  tne 
interrupt  vector  address  are  selectable  with  switches  on  the  quad 
board  and  are  currently  set  at  octal  767700  and  0  respectively. 
(Although  interrupts  are  not  used  now,  they  are  available  for  future 
use.)  The  defined  addresses  and  their  descriptions  are  as  follows: 


5ISML  APgRESS  0E3SBIEIIQH 

CSR  Base  Controi/status  word. 

Bit  0:RL?.ST  -  A  zero  in  this  bit 
stops  HDG  controller  and  sets 
program  to  address  0.  A  one  enables 

HDG  operation.  Read/write. 

Bit  1:Fage  -  Selects  program  from 
either  page  0  or  page  1  of  ROM. 
Read/write . 

Bit  2: Not  used. 

Bit  3:STEXP/  -  Starts  exposure 
sequence  when  raised  to  a  one. 

Should  be  held  at  zero  until  all 
descriptors  have  been  passed  from 
LSI  11,  then  raised  for  at  least  1 
instruction  cycle.  Must  be  returned 
to  zero  before  next  DMA  transfer  is 
complete.  Read/write. 

Bit4:Not  used. 

Bit  5:ENAST  -  Enable  interrupt  A. 
Read/write. 


'  Bit  6 :ENBST  -  Enable  interrupt  3. 
Read/write . 


Bit  7:RGSTA  -  Interrupt  request  A. 

Not  currently  used.  Read  only. 

Bit  8:ECC  -  End  of  cell  bit  indicates 
when  one  for  more  than  10  microsec¬ 
onds  that  a  complete  cell  has  been 
written.  May  be  one  for  less  than  10 
microseconds  at  end  of  individual 
scan  lines.  Read  only. 

Bits  9  through  i4:Not  used. 

Bit  15:RQSTB  -  Interrupt  request  B. 
Not  currently  used.  Read  only. 


XCFF  Base+2  X  offset  in  pixels  in  bits 

0-11  in  two’s  complement  format. 
Write  only 


I OFF  Base+4  Y  offset  in  scan  lines  in 

bits  0-9  in  positive  true  format. 
Write  only 

YMAG  Base+lQ  Y  magnification  factor  in 

bits  C-7.  Value  of  128  (decimal) 
represents  nominal  magnification 
factor.  Values  of  255  and  0  repre¬ 
sent  approximate  changes  of  +/-  10 
percent.  Write  only 


L.RPT  Base*  12  Number  of  scan  lines  per 

bar  group. LSPT  equals  two’s  com¬ 
plement  of  (scan  lines  ♦  1).  Bits  0-9 
only.  Write  only. 

XMAG  Base+14  X  magnification  factor  in 

bits  0-7.  Write  only. 

The  DMA  interface  consists  of  a  DEC  GRV118  interface  board  in  the 
computer  and  complementary  circuitry  on  the  HDG  board.  As  with  the  PIO 
interface,  the-two-  are- connected  by  two  20  foot  ribbon  cables. 

The  DMA  base  address  can  be  set  with  switches  on  the  DRV  1 1 B  board  and 
is  currently  set  at  octal  772410.  The  registers  in  the  DRV  1 1 B  of 
interest  here  are  the  word  count  register  (WCR)  at  the  base  address, 
the  bus  address  register  (BAR)  at  the  base  address  plus  2,  and  the 
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control/status  register  (CSR)  at  the  base  address  plus  Details  of 
DRV  1 1 B  operation  can  be  found  in  the  appropriate  DEC  documentation; 
knowledge  of  this  information  is  assumed  in  the  following  discussion. 

The  LSI  11  computer  signals  the  HDG  when  it  ready  to  make  a  DMA 

transfer  by  dropping  the  READY  line.  The  HDG,  when  it  reaches  the 
appropriate  point  in  its  cycle,  examines  the  READY  line  and,  as  soon 
as  it  is  found  low,  starts  the  DMA  transfer.  The  transfer  proceeds  in 
the  single  cycle  mode  as  fast  as  the  computer  will  permit. There  is  no 
further  intervention  by  either  the  computer  or  the  HDG  until  the  WCR 
overflows,  thus  raising  the  READY  line  and  terminating  the  operation. 
With  the  current  version  of  the  HDG  program,  the  READY  line  is 
monitored  by  the  HDG  and  other  tasks  are  suspended  until  it  rises.  A 
transfer  of  1024  words  (  the  size  of  the  HDG  memory)  takes  about  3.5 
milliseconds 

The  data  are  transferred  to  the  HDG  in  a  word  format  and  input 
directly  to  the  HDG  memory.  The  HDG  interface  generates  memory 
address  increment  and  memory  write  signals  automatically. 

Although  only  one  way  data  transfer  is  required  (computer  to  HDG), 
there  is  provision  for  transferring  information  from  the  HDG  memory 
back  to  the  computer  for  diagnostic  purposes.  This  is  done  by  raising 
the  FNCT  1  bit  (bit  1)  in  the  DRV  1 1B  CSR. 

Circuitry  is  included  for  testing  single  word  parity  during  a  computer 
to  HDG  transfer.  The  parity  bit  is  bit  13  of  each  word  and  may  be 
either  even  or  odd  depending  on  the  condition  of  jumpers  on  the  HDG 
circuit  board.  The  jumpers  are  currently  set  for  even  parity.  The 
parity  test  is  selected  by  raising  the  FNCT  2  bit  (bit  2)  in  the 
DRV  1 1 B  CSR.  Upon  detection  of  a  parity  error,  the  HDG  will  raise  the 
ATTN  signal  (bit  13)  which  will  in  turn  set  the  ERROR  bit  (bit  15), 
terminate  the  transfer,  and  generate  an  interrupt  if  the  interrupt 
enable  bit  (bit  6)  is  set. 

The  remaining  function  and  status  lines  available  on  the  DRV11B  are 
not  used. 


3.2  Memory 

The  memory  is  arranged  in  a  1024  x  16  bit  configuration  and  includes  a 
memory  address  counter  and  an  address  return  register.  The  address 
counter  can  be  reset, preloaded  from  the  return  register,  or 
incremented  by  a  HDG  controller  instruction.  The  counter  is 
incremented  automatically  during  a  DMA  operation  after  each  word 
transfer.  The  return  register  can  be  loaded  with  the  contents  of  the 
address  counter  with  a  controller  instruction.  Write  signals  to  the 
memory  are  generated  by  the  interface  circuitry  during  a  DMA 
operation. 

Bar  widths  are  stored  sequentially  starting  at  location  0.  Bits  0-11 
of  each  word  represent  bar  width  in  a  two's  complement  format.  Bit  15 
is  a  1  if  the  corresponding  bar  is  to  be  unblanked.  The  cumulative 
value  of  bar  widths  in  any  group  must  total  at  least  2400  pixels  less 
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V 


che  value  of 


tne  X  offset.  The  cumulative  value  of  ail  the  bar  widths 
except  the  last  in  a  group  must  be  less  than  2000  pixels.  Meeting 
these  criteria  ensures  that  the  last  descriptor  in  a  group  will  be 
active  at  the  end  of  eacn  scan  line.  Failure  to  meet  these  criteria 
will  result  in  misregistration  between  bar  groups. 

The  last  word  in  a  block  of  ceil  descriptors,  following  the  last  oar 
widtn  in  the  last  bar  group,  signifies  the  end  of  ceil  with  a  1  in  bit 
14. 

The  remaining  bits  (12,13)  are  not  used. 


3.3  X  Deflection 

The  X  deflection  signal  is  generated  by  an  integrating  sawtooth 
generator  synchronized  with  a  counter.  The  counter  circuit  has  2 
counter  chains  which  are  driven  by  a  stable  10  MHz  clock  and  count  up 
to  0  from  preset  values.  Additional  circuitry  generates  2  signals, 
called  S Vi  and  XVAL,  wnich  are  used  tc  synchronize  the  other  parts  of 
the  HDG.  The  SW  signal  controls  solid  state  switches  in  the  feedback 
path  of  an  integrating  amplifier.  When  SW  is  1 ,  the  switches  are  open 
and  the  integrator  generates  a  ramp  waveform.  When  the  switches  are 
closed,  the  waveform  retraces  rapidly.  SW  is  high  for  36  microseconds 
and  low  for  16  microseconds.  These  values  can  be  changed  in  100 
nanosecond  increments  by  changing  jumpers  on  tae  HDG  circuit  board. 

The  XV  aL  signal  rises  5  microseconds  after  the  SW  signal  and  falls 
simultaneously.  XVAL  defines  tae  useful  sweep  time.  The  6 
microsecond  hiatus  between  SW  ar.c  XVAL  provides  time  for  the  sweep 
waveform  to  become  stable  and  linear. 

Because  deflector,  repeatability  is  important,  the  X  sweep  runs 
continuously  and  covers  the  full  2400  pixel  format.  Offset  of  the 
origin  in  the  X  direction  is  achieves  by  appropriate  blanking  at  the 
beginning  and  end  of  each  line  rather  than  by  moving  the  entire 
pattern. 

The  reference  signal  used  tc  drive  the  integrator  is  provided  by  a 
precision  reference  voltage  conditioned  by  the  XMAG  word  via  an  S  bit 
D/A  converter.  The  result  is  a  reference  signal  which  can  be  varied 
by  the  XMAG  word  by  +/-  10  percent  in  steps  of  1  part  in  256. 

Trim  controls  have  been  provided  to  vary  the  nominal  scale  factor,  the 
gain  of  the  XMAG  word,  the  nominal  X  center  value,  and  the  value  of 
ihe  reference  voltage. 


3.4  Y  Deflection 


The  Y  deflection  circuitry  includes  a  10  bit  line  counter,  a 
D/A  converter,  magnification  control  circuitry,  and  an 
amplifier. 


10  bit 
output 


The  counter  can  be  reset,  incremented,  or  loaded  with  the  Y 


offset 


D-8 


word  from  the  LSI  11  under  microprogram  control.  With  a  range  of  10 
Pits,  up  to  1024  raster  lines  can  be  addressed.  The  pattern  is  offset 
to  accomodate  positioning  tolerances  by  preloading  the  counter  with 
the  Y  offset  word  (YOFF).  Subsequent  lines  are  addressed  by 
incrementing  the  counter  one  line  at  a  time. 

The  D/A  converter  is  a  12  bxt  unit  with  the  2  least  significant  bits 
tied  to  zeroes.  As  with  the  X  sweep,  the  reference  drive  is  derived 
from  a  precision  voltage  source  conditioned  by  the  YMAG  word.  Full 
scale  variation  of  the  YMAG  word  will  cause  a  +/-  10  percent  change  in 
the  reference  signal  in  steps  of  1  part  in  256. 

The  ouput  amplifier  is  the  same  as  in  the  X  axis.  It  is  capable  cf 
driving  a  50  ohm  load  at  */-  5  volts  and  has  a  bandwidth  in  excess  of 
50  MHz. 

Trim  controls  have  been  provided  to  vary  the  nominal  pattern  size, 
position  offset,  magnification  factor,  and  reference  voltage. 


3.5  Unblanking 

The  CRT  image  is  generated  by  converting  the  bar  width  descriptors  to 
time  delays  in.  synchronism  with  the  X  sweep  and  using  the  resulting 
signal  to  turn  the  CRT  beam  on  and  off. 

The  bar  width  conversion  circuitry  includes  2  counter  chains  for  the 
coarse  timing  and  a  variable  delay  circuit  for  the  fine  timing.  The 
bar  width  descriptors,  which  are  stored  sequentially  in  the  memory, 
are  12  bit  numbers  in  two’s  complement  notation  specifying  the  numDer 
of  pixels  in  each  bar.  The  2  least  significant  bits  control  a  delay 
circuit  which,  at  the  appropriate  time,  generates  a  delay  of  0,  12.5, 
25,  or  37.5  nanoseconds  corresponding  to  0  through  3  pixels.  The 
remaining  10  bits  are  used  to  preset  the  counters. 

The  counters  are  used  alternately;  while  one  is  counting,  the  other  ia 
being  loaded  with  tit  next  bar  descriptor.  The  active  counter  is 
selected  by  an  EVEN/ODu  flip  flop  that  is  toggled  at  the  end  of  each 
count  sequence,  thereby  transferring  control  to  the  opposite  counter. 
Both  counters  are  enabled  or  disabled  simultaneously  by  the  CNTEN  flip 
flop  which  is  set  or  cleared  by  the  controller.  The  counters  are 
clocked  by  a  20  MHz  clock  and  advance  from  their  preset  values  till 
they  overflow  and  a  carry  signal  is  generated.  The  carry  signal 
resets  the  counter,  triggers  the  delay  circuit,  toggles  the  EVEN/ODD 
flop,  and  sets  a  status  flip  flop  monitored  by  the  controller.  The 
opposite  counter  is  enabled  immediately  and  proceeds  to  count  on  the 
next  rising  edge  of  the  20  MHz  clock. 

The  2  least  significant  bits  have,  prior  to  the  advent  of  the  carry 
signal,  been  buffered  in  latches  and  decoded  to  select  1  of  4  possible 
delay  paths.  The  carry  signal  is  propagated  through  the  selected  path 
and,  at  the  other  end,  clocks,  into  a  flip  flop,  bit  15  of  the  bar 
descriptor,  corresponding  to  the  value  of  the  next  bar,  i.e.  1  or  0 
indicating  blanked  or  unblanked. 

The  unblanking  signal  at  the  output  of  the  flop  is  buffered  and 
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converted  to  a  voltage  swing  of  0  to  -2-5  volts  with  an  iipedance  of 
50  ohms  to  drive  the  Z  input  of  the  Infodex  display. 


When  the  controller  detects  a  carry  signal  indicating  counter 
overflow,  it  executes  the  steps  necessary  to  prepare  the  disablea 
counter  for  the  next  sequence.  This  includes  incrementing  the  memory 
and  loading  the  counter  and  delay  Puffer.  These  activities  require 
about  400  nanoseconds  (4  instruction  cycles)  so  bar  widens  smaller 

tr.an  32  pixels  (12.5  r.3  per  pixel)  cannot  be  accommodated. 


When  one  counter  overflows,  the  other  starts  counting  on  the  next  20 
KHz  clock  cycle  irrespective  of  the  value  of  the  delay.  To  prevent 
truncation  error,  which  will  accumulate  if  not  corrected,  a  modified 
bar  width  is  generated  by  subtracting  the  time  delay  associated  with 
the  previous  bar  from  the  current  bar  width.  It  is  this  modified 
winch,  rather  than  the  value  directly  from  memory,  which  is  directed 

co  the  counter  and  delay  circuits. 


5.6  Line  counter 

In  a  given  hologram,  all  bars  are  a  predetermined  size  in  the  Y 
dimension.  This  size  is  specified  by  the  Line  Repeat  word  (LRPT) 
transmitted  from  the  LSI  1 1  via  the  PIG  interface  and  maintained  in  a 
10  bit  counter.  The  line  repeat  counter  can  be  reset,  incremented,  or 
loaded  with  the  LRPT  word  by  the  controller.  The  carry  bit  of  the 
counter  can  be  monitored  by  the  controller  to  indicate  overflew. 
Because  the  counter  cannot  be  decremented.  LRPT  must  be  in  a  two's 
complement  format. 


3.7  Controller 

The  hDG  controller  operates  on  a  32  bit  microword  in  accordance  with  a 
PRCM  stored  program  at  a  cycle  rate  of  '.0  MHz.  It  can  issue  control 

signals,  monitor  the  status  of  input  signals,  and  branch  in  response 
to  programmed  conditions. 

Two  types  of  branch  are  permitted:  a  tabular  branch  and  an 
unconditional  branch.  In  the  unconditional  branch  mode,  when  a 
selected  status  signal  becomes  true,  the  program  branches  to  a 
previously  specified  location.  In  the  tabular  branch  mode,  the  bit 
pattern  formed  by  the  condition  of  a  group  of  status  signals  is  OR'd 
with  a  base  address  to  yield  a  branch  location. 

The  microword  is  divided  into  4  groups  of  3  bits  each.  Bits  0  through 
7  specify  the  address  of  the  next  instruction  or  the  base  address  of  a 
tabular  branch  instruction. 

Bits  8  through  15  are  mask  bits  used  to  select  status  signals  for 
examination.  The  lower  5  bits  of  the  group  are  AND'd  with  5  status 
signals  to  form  a  tabular  branch  word. 


BITS 

0 

SISM1 

3P  V  « 

DESCRIPTION 

Carry  from  counter  * .  Must  be 

cleared  by 

- 

CRY2 

CLRCRY . 

Carry  from  counter  2.  Must  be 

cleared  zy 

- 

V  T 

CLRCRY. 

C^rry  c  c  ur  t  ^ r%  ^  v *  —  & 

cl1?2  r?- 

3 

CRY 

CLRCRY. 

Carry  from  any  of  above.  Must 

be  cleared 

ii 

READY 

by  CLRCRY. 

Goes  high  when  DMA  READY  line 

goes  low. 

Table  1 

-  Tabular 

3ranch  Cedes 

Thu  upper  3  bi 

ts  are  decoded  to  select  1  of  7  status  signals 

(or  none 

of  the  7)  to  be  monitored 

for  an  unconditional  branch.  Both 

types  of 

branching  can 

be  used  in 

a  single  instruction.  If  the  unconditional 

branch  status 

is  false. 

the  tabular  branch  is  executed. 

If  the 

unconditional 

branch  stab 

us  becomes  true,  the  unconditional 

Branch  is 

executed . 

CODE 

SIGNAL 

DESCRIPTION 

20 

RUN 

false  when  DMA  is  active 

u0 

LNCRY/ 

True  when  line  repeat  counter 

overflows 

60 

XVAL 

false  during  X  sweep 

60 

XVAL/ 

True  during  X  sweep 

AO 

YCRY/ 

True  when  Y  sweep  counter 

overflows 

CO 

STEXP/ 

True  when  bit  3  of  CSP. 

is  true 

E0 

EOC/ 

True  when  End  of  Cell  line 

(bit  14  of  bar  width  descriptor) 
is  true 


Table  2  -  Unconditional  Branch  Codes 


Sits  16  through  23  define  function  codes  which  are  decoded  to  generate 
control  signals.  The  upper  and  lower  nibbles  are  decoded  separately 
into  15  lines  each  plus  NOP 'a.  By  specifying  functions  in  both 
nibbles,  2  control  signals  can  be  generated  in  a  single  100  nanosecond 
instruction  cycle. 


u-1 1 


CQD1 


hMEMSIia 


i 


0 

NOP 

No  operation 

1 

MR 

Maser  reset 

2 

LDCNTC 

Load  ode  counter 

3 

LDCNTS 

Load  even  counter 

4 

LDDLY 

Load  delay  buffer 

5 

LDLC 

Load  line  repeat  counter 

6 

INCL 

Increment  line  repeat  cou; 

7 

RSTL 

Reset  line  repeat  counter 

3 

BMSET 

Enable  beam  on 

9 

STEVEN 

Set  EVEN/ODD  flag  to  even 

A 

DMAST 

Start  a  DMA  transfer 

Ma 

1 

0 

NOP 

No  operation 

i 

MARST 

Memory  address  reset 

2 

MAENB 

Store  memory  address 

3 

LDM 

Recall  memory  address 

4 

INCM 

Increment  memory  address 

5 

CNTON 

Enable  counters 

6 

CNTOFF 

Disable  counters 

7 

Not  used 

8 

LOAD 

Load  branch  address 

9 

YRST 

Reset  Y  counter 

A 

YLD 

Load  Y  counter  with  YOFF 

3 

YC'JT 

Increment  Y  counter 

C 

LDXOFF 

Load  X OFF  into  counters 

D 

CLRCRY 

Clear  the  carry  flags 

Table 

3  -  Function 

Codes 

Bits  24  through  31  specify  the  location  of  the  branch  address  in  an 

unconditional  branch  instruction.  The  branch  address  can  be  issued  in 
any  instruction  before  or  including  the  branch  instruction. 

The  program  memory  has  a  capacity  of  512  words  divided  evenly  into  2 
pages.  The  page  is  selected  by  bit  1  of  the  HDG  CSR  word.  Controller 
programs  cannot  operate  over  the  page  boundary. 

The  controller  is  enabled  by  RLRST  (bit  0  of  the  HDG  CSR  word).  When 
RLRST  is  0,  the  controller  is  disabled.  When  RLRST  is  driven  to  a  1 , 
the  controller  starts  operating  at  location  0  of  the  selected  page. 
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4.0  MICROCODE 


Following  are  flowcharts  and  coding  sheets  describing  the  HDG 
microcode.  Each  block  on  the  flowchart  corresponds  to  a  single 
instruction.  The  numbers  to  the  left  of  the  blocks  refer  to 
instruction  addresses. 

The  program  has  several  nested  loops  which  are  executed  with  varying 
frequency.  Instruction  O'  is  encountered  only  a”  initial  turn  on. 
Instructions  02  through  OB  (between  entry  points  J  and  i)  are  executed 
once  per  cell.  The  loop  starting  at  0C  occurs  once  every  n  lines 
where  n  corresponds  to  the  specified  number  of  line  repeats.  The  loop 
starting  at  0E  is  executed  every  line.  Most  of  the  action  takes  place 
between  entry  points  2  and  3;  this  loop  is  traversed  once  every  two 
bar  widths  which,  at  maximum  -esolution,  results  in  a  loop  cycle  time 
of  about  I  microsecond. 

After  a  master  reset  and  a  memory  address  reset,  the  controller  waits 
for  the  DMA  READY  line  to  fail  indicating  that  the  LSI  1 1  is  ready  for 
a  transfer.  During  the  'wait,  the  3MRST  command  is  continually  issued. 
When  the  READY  line  fails,  the  counters  are  turned  off,  the  memory 
address  is  reset  again,  the  transfer  is  started,  and  the  branch 
address  for  the  next  instruction  is  loaded.  The  controller  then 
monitors  the  RUN  line  to  determine  when  the  transfer  is  complete. 
After  the  transfer,  the  program  loops  until  the  STEXP/  line  falls 
indicating  the  LSI  is  ready  for  the  exposure  sequence. 
Instructions  09  and  0A  set  the  EVEN/CDD  flag  to  even  and  load  the  X 
and  Y  offsets.  Instructions  33  through  3C  synchronize  the  program 
with  the  X  sweep.  The  following  instruction,  CB,  is  executeed 
immediately  after  the  start  of  the  retrace. 

The  loop  starting  at  0C  is  executed  for  the  specified  number  of  line 
repeats.  At  the  3tart  of  the  loop,,  the  memory  address  is  saved  so  it 
may  be  recalled  at  the  start  of  each  line  within  the  group,  and  the 
line  repeat  counter  is  reloaded  with  the  line  repeat  count.  The  loop 

starting  at  0E  is  executed  on  each  scan  line.  The  timing  is  such  that 
this  instruction  occurs  during  the  beginning  of  the  retrace  period. 
So.  after  loading  the  X  offset,  the  program  loops  until  the  start  of 
the  sweep. 

The  counters  are  immediately  enabled  and  the  unblanking  loop, 
comprising  instructions  ’A  through  22,  is  entered.  This  is  a 

continuous  loop  and  is  exited  only  when  XVAL  falls  indicating  the  end 
of  the  forward  sweep  period  and  the  start  of  the  retrace.  This  signal 
is  examined  in  instructions  1C  and  20  independently  of  other  loop 
operations.  The  remaining  instructions  set  up  the  counters  and  delay 
circuit  to  measure  out  the  bar  widths  and  shuttle  operation  between 
the  two  counters.  At  each  bar,  the  program  loops  until  the 
appropriate  counter  overflows,  then  clears  the  carry  flag,  loads  the 
deiay  for  the  current  bar,  incrments  the  memory  address,  and  loads  the 
opposite  counter.  The  sequence  is  repeated  as  many  times  as  necessary 
until  the  scan  Line  13  complete. 

At  the  end  of  each  bar,  the  counters  are  turned  off  and  the  line 
repeat  and  Y  position  counters  are  incremented.  The  LNCRY/  line  is 
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tested  to  determine  if  the  line  repeat  counter  has  overflowed.  If 

not,  the  previously  stored  memory  address  is  reloaded  into  the  memory 
address  counter  and  operation  proceeds  from  instruction  OE. 

If  the  line  repeat  counter  has  overflowed,  bit  14  of  the  next  word  in 

memory  is  tested  to  determine  if  the  end  of  the  cell  has  been  reached. 
If  not,  the  program  branches  to  OC  to  generate  the  next  group  of 

lines.  If  it  has,  control  returns  to  instruction  02  to  await  the  next 
DMA  transfer. 
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MISSION 
of 

Rome  Air  Development  Center 

RAVC  pZa.ni>  and  execute*  research,  de.veZopme.nt,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control 
Cormunccationi  and  Intelligence  (Ch)  activities.  Technical 
and  engA.neeAi.ng  support  within  areas  of  technical  competence 
AA  provA.ded  to  E SO  Program  Offices  IPOs)  and  other  BSD 
elements.  The  pAlncipal  technical  mission  areas  ane 
commu.nicatA.ons ,  electromagnetic  guidance  and  control,  sur- 
velltance  of  ground  and  aerospace  objects,  intelligence  data 
cou.eatA.on  and  handling,  information  system  technology, 
A.onospherA.c  propagatcon,  solid  state  sciences,  microwave 
physios  and  electronic  reliability,  maintainability  and 
compatibility. 
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